Scaling of the Hysteresis Loop in Two-dimensional Solidification 
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The first order phase transitions between a two-dimensional (2d) gas and the 2d solid of the first 
monolayer have been studied for the noble gases Ar, Kr and Xe on a NaCl(fOO) surface in quasi- 
equilibrium with the three-dimensional (3d) gas phase. Using linear temperature ramps, we show 
that the widths of the hysteresis loops of these transitions as a function of the heating rate, r, scales 
with a power law oc r a with a between 0.4 and 0.5 depending on the system. The hysteresis loops 
for different heating rates are similar. The island area of the condensed layer was found to grow 
initially with a t 4 time dependence. These results are in agreement with theory, which predicts 
q = 0.5 and hysteresis loop similarity. 
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Adsorption isotherm measurements are long known to 
give detailed information about adsorbated-substrate in- 
teractions as well as about lateral interactions, if com- 
bined with structural information. They were first intro- 
duced into surface studies by Langmuir |l]] . In most cases, 
adsorption is a first order phase transition || with the 
heat of adsorption as the latent heat, which is expected 
to show the typical hysteresis behavior due to overheat- 
ing and undercooling the adsorbed layer. In this Letter 
we want to concentrate on a specific aspect of hysteresis, 
scaling of the energy turn-over as a function of frequency 
and amplitude, as measured by the hysteresis loop area. 

Scaling of the hysteresis loop area by cycling through 
a first order phase transition has been predicted long ago 
and has been studied in detail in ferromagnetic systems. 
Already in early pioneering work on a 3D magnet, which 
dates back to the last century |3| , first indications of scal- 
ing behavior of the hysteresis loop area have been found. 
Recently, experimental investigations of hysteresis loop 
scaling in ultrathin magnetic films have been performed 
|, ||, ||. However, although power law behavior is 
found in all studies [Q, ||, g, the effective exponents 
obtained differ by typically one order of magnitude be- 
tween theoretical models Q and measurements of energy 
loss scaling in ultrathin magnetic films (see discussion in 
Ref. ||). No obvious explanation for this large discrep- 
ancy has come up yet. 

Universality mean that scaling of the hysteresis loop 
area to be observable not only in ferromagnetic systems, 
but also, e.g., in adsorbed layers, for which, according to 
our knowledge, this property has not been studied in any 
detail yet. This class of systems is characterized not only 
by completely different parameters, but also by narrow 
domain walls of just a few lattice constants, in contrast 
to magnetic systems. As a result nucleation barriers in 
adsorbed layers can be much smaller than in magnetic 
films. We have recently predicted scaling behavior of 
the adsorption hysteresis loop width, fi^, and of the hys- 



teresis loop area A || as a function of amplitude and 
cycle frequency for the situation close to thermodynamic 
equilibrium. A more detailed analysis of the hysteresis 
behavior for the case of magnetic film jl0| predicts 
also similarity of different hysteresis loops. 

For our test experiments of scaling of hysteresis loops 
in adsorption of two-dimensional layers, we used the 
phase transition to condensation of a 2d solid of the no- 
ble gases Xe, Kr and Ar on single crystalline thin films of 
NaCl(lOO) in quasi-equilibrium with the 3d gas phase at 
room temperature. Due to the van-der-Waals attraction 
between the noble gas atoms, this transition is strongly 
first order, and can be easily measured, as shown be- 
low. The deviation of the chemical potential from the 
equilibrium value at the phase transition, fi, in the case 
of adsorption plays the role of the magnetic field, which 
was periodically varied around the equilibrium value with 
amplitude /xo- 

In our model, the process of adsorption is divided into 
two stages: nucleation of islands and growth of these is- 
lands. The situation when the nucleation time t„ -C t, 
the cycle time, has been considered in Ref. ||. We call 
this situation growth controlled hysteresis. For heteroge- 
neous nucleation this condition can be easily fulfilled. 

With the additional assumption that an island with 
area A oc L 2 grows with a boundary velocity that 
depends on /i like v — Tfi, we obtain for the island of 
diameter L 



L(t)-L = 2CrTfi (^P-) 2 



CT 



(1) 



where Lq is of the order of the size of a stable nucleation 
center, ji(t) — rt with linear heating rate r (adapted ac- 
cording to the experiments carried out), and C a constant 
of order unity. We have explicitly used that (t„/t) 2 <C 1. 
Eq.|l| shows that the system behavior is controlled by the 
characteristic length L T = tT/j,q = A[± 2 Y/r. With in- 
crease of L T the area of the adsorption-desorption cycle 
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loop will decrease. If Ln is the mean distance between 
nucleation centers, then the hysteresis parameters are de- 
fined by the ratio L N /L T §, f§. 

From Eqjl] the scaling behavior of the hysteresis loop 
width fih and of the area of the adsorption-desorption 
cycle loop A as functions of r and no can be evaluated 
in the limit Ln *C L Tl i.e. when Ln determines the 
maximum island size. At fi = fih/2 the typical island size 
L w Ln. The condition Lpj <C L T requires that \ih *C Mo 
and Lff (rr//io)(/i/ l ) 2 . If we assume that Lm does not 
depend significantly on r and /io, we immediately obtain 
that /ih scales as [ih oc yjUo/r oc y 7 ?. Under the same 
conditions of Ln independent (or weakly dependent) on 
Ho and r, the hysteresis loops are similar and scale in 
the same way as fjt,^. This has been predicted in 0. At 
sufficiently low frequencies, the height of hysteresis loop 
does not depend on frequency. 

This part of the frequency dependence is tested in our 
experiments, in which we linearly ramp the surface tem- 
perature up and down, keeping the ambient gas pressure 
constant. In addition to scaling of the hysteresis loop 
area, the time dependence of island sizes can be tested 
by experiment, which, according to eq. [I] should grow 
as L(t) oc rt 2 . As long as the islands don't touch each 
other, the coverage 9 depends on L as oc L 2 oc r 2 i 4 . 
In a diffraction experiment, the integrated intensity of a 
superstructure spot should be proportional to the cover- 
age so that this equation can be tested easily measuring 
the initial time dependence of the integrated diffracted 
intensities. 

Our experiments were carried out in a UHV-chamber 
equipped with high resolution LEED supplemented by a 
high resolution optical detector jl]J , Auger spectrometer 
and a quadrupole mass spectrometer, at a base pressure 
of 2 x 10~ n mbar. Peak intensities of LEED diffraction 
spots were both measured with a Faraday cup and with 
the optical detector, while integrated intensities were ob- 
tained from images taken from the optical detector using 
a high resolution slow scan CCD camera. The NaCl(lOO) 
surfaces were prepared in situ as 3 double layers thick epi- 
taxial films by evaporation of NaCl onto a Ge(100) sub- 
strate at a surface temperature of 200 K and subsequent 
annealing to 550 K. This procedure produces single crys- 
talline NaCl layers that are in registry with the Ge sub- 
strate on the flat terraces and overgrow monoatomic steps 
of the Ge substrate in a carpet-like mode JlJ] . Therefore 
the quality of the films is determined by the step den- 
sity of the substrate, which had typical terrace length 
of 400 A. The Ge substrates (size 10 x 10 x 1mm) were 
oriented with a precision of 0.1° using a Laue camera 
and polished with diamond pastes and a final chcmo- 
mechanical treatment. These samples were mounted onto 
a He-cooled flow-cryostat. They were heated either by 
electron bombardment or radiation. For the measure- 
ments only radiation was used. Temperature was mea- 
sured using a Chromel-Alumel thermocouple attached to 
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FIG. 1. Hysteresis loops of the 2d gas-solid phase transi- 
tion of Xe/NaCl(100) at a Xe pressure of 1(T 7 mbar. The 
peak intensity of the of a first order superstructure spot was 
monitored at the heating rates indicated. 

the base plate of the mounting assembly, which was in 
good thermal contact with the Ge substrate. A comput- 
erized temperature controller with a nominal resolution 
of 0.001 K was used for linear temperature ramps up and 
down. High purity gases were dosed directly onto the 
sample through a ring-shaped slit centered around the 
surface normal, which could be cooled to 100 K with |N2. 
Measurements carried out both with gas at room tem- 
perature and at 100 K showed that the gas temperature 
had a negligible influence on the results presented be- 
low. Pressures given below are pressures obtained from a 
calibration with gas exposures from background gas. Ho- 
mogeneity of the pressure on the surface was also tested 
and found to vary negligibly over the surface area seen 
by the LEED beam. 

Xe on NaCl(100) condensed in the first monolayer 
forms a quasi-hexagonal incommensurate structure |l3| ] 
so that superstructure diffraction spots were directly used 
to measure hysteresis both with integrated and peak in- 
tensities. While the first monolayer of Ar forms an or- 
dered (1 x 1) structure, a diffuse (2 x 1) structure with 
glide plane symmetry is seen for Kr. Only peak intensi- 
ties of integral order beams have been evaluated for the 
latter two systems. Data of equivalent beams have been 
averaged where available. 

Typical data of the hysteresis during condensation of 
the quasi-hexagonal Xe layer are shown in Fig. [I]. The 
linear heating rate was varied by two orders of magnitude 
between 0.001 and 0.1 K/s. This rate was limited at small 
rates by the resolution of temperature measurement, at 
high rates by the onset of intermixing of second layer con- 
densation. Measurements were taken for gas pressures of 
10~ 6 and 10~ 7 mbar. Not surprisingly, slopes during ad- 
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FIG. 2. Hysteresis loops of the 2d gas-solid phase transition 
of Ar/NaCl(100). The change of peak intensity of an first 
order integral spot was monitored for the 2d-(l x 1) Ar solid 
formed. 

sorption and desorption are not symmetrical due to the 
exponential dependence of the desorption rate on sur- 
face temperature. At the highest rate the form during 
condensation changes due to the mentioned onset of in- 
termixing with second layer adsorption, leading also to a 
decrease of intensity with decreasing temperature. Inside 
these limits of heating rate the intensity at saturation did 
not depend on the heating rate. 

Similar measurements were taken measuring the con- 
densation induced changes of peak intensities of integral 
order beams for Ar and Kr on NaCl, which, depending on 
electron energy of the LEED experiment and on diffrac- 
tion order, can be both positive or negative. An exam- 
ple for Ar induced intensity changes is shown in Fig. |^. 
Please note that the dependence on noble gas concentra- 
tion need not be linear in these cases. 

For all three systems, we evaluated the dependence of 
the widths of the hysteresis loops at half maximum in- 
tensity on the heating rate r. The results are shown as 
a log-log plot in Fig. ^. Power laws were obtained for 
all three systems. The effective exponents a obtained 
from these plots are close to 0.4 for Xe and Ar conden- 
sation (within the statistical uncertainty of about 10%), 
whereas the average value for Kr of a = 0.50 actually cor- 
responds exactly to the value expected from the simple 
model of growth controlled hysteresis already mentioned. 
Also the experimental results of Xe and Ar are sufhcently 
close to this value so that this model seems to describe 
the essential physics correctly. 

A further test of this model can be carried out by ana- 
lyzing the time dependence of island growth in the initial 
stages of growth. For this purpose we plotted the inte- 
grated intensities of the Xe systems as a function of time 
during condensation again on log- log scale (see Fig El) . As 



FIG. 3. Log-log plots of the width of hysteresis loops, 
AT, as a function of the heating rate r in presence of a 
3d gas pressure. Xe/NaCl(100): p = 1 x l(T 7 mbar (O), 
p = 1 X lCrVbar (□). Full symbols: Kr/NaCl(100) at 
p = 1 X 10~ 7 mbar (three different orders of diffraction), o 
and A: Ar/Nacl(100) at p = 1 x 10~ 7 mbar measured for the 
(10) and (11) beams, respectively. 

mentioned, the model predicts the integrated intensity to 
increase oc r 2 i 4 , i.e. cx (AT) 4 /r 2 . Therefore, the data in 
this figure are plotted once as a function of AT, and a 
second time scaled by r 2 . The integrated intensities fol- 
low closely the t 4 dependence predicted by our simple 
model, but deviates from it when coverage gets closer to 
saturation. Of course, this plot is very sensitive to the 
choice of time zero, for which we chose the condensation 
temperature at equilibrium, as estimated from the cen- 
ter of the hysteresis curves. This is the earliest possible 
time. It is fully consistent with the assumptions of the 
model used, which assumes a small nucleation time r„. 
The intensity zero was taken as the bottom of the hys- 
teresis curves without further adjustments. This result 
therefore seems to nicely corroborate the model assump- 
tion of growth controlled hysteresis, which is not limited 
by diffusion on the surface. 

These results turned out to be insensitive to additional 
production of anion vacancies on the NAC1 substrate, 
which is caused by the measuring electron beam. They 
act as additional nucleation centers. Though they re- 
duce the maximum size of islands, they obviously do 
not change the growth modes, in agreement with the ex- 
pectations from our model. For larger islands, diffusion 
to the island boundaries of course can no longer be ne- 
glected resulting in deviations from our model. Diffusion 
is particularly important in the second layer on already 
condensed islands, since there the sticking coefficient is 
much larger than on the bare surface jL2j. This process 
increases the speed of growth, especially at the initial 
stages. Although in a real experiment mass transport to 
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FIG. 4. Test of scaling with time of integral intensities of a 
first order superstructure spot of Xe/NaCl(100), The curves 
close to the straight line have been divided by r 2 . 
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the growing boundary by diffusion is always present, we 
therefore expect that our main result -dynamical scaling- 
remains valid for conditions close to equilibrium. 

An explicit test of similarity of the hysteresis curves 
was again carried out for Xe/NaCl(100). The rescaled 
curves (after centering) are shown in Fig. || for one or- 
der of magnitude changes in the heating rate. While the 
evaporation data fit perfectly to a common line, there is 
more scatter in the data during condensation but no gen- 
eral trend for the small deviations was found. Therefore, 
also similarity seems to be fulfilled by these data. 

To conclude, we have introduced a new method of 
studing adsorption phenomena, observed scaling behav- 
ior for adsorption, found corresponding critical exponents 
and proved similarity of hysteresis loops. Our studies of 
scaling of the hysteresis loop areas in adsorption gave re- 



sults in almost quantitative agreement with theoretical 
expectations using a simple model. This situation is in 
sharp contrast with the situation in magnetic films. Pos- 
sible reasons are the small width of the interface between 
islands (of the order of one lattice constant) and/or a low 
interface energy. The latter would cause small nucleation 
barriers. It is consistent with the dominant mechanism 
of growth controlled hysteresis found in our experiments. 
For studies of hysteresis phenomena, we have shown that 
adsorbed films are complementary to ultrathin magnetic 
films, since they open a completely different parameter 
space. The method we have used can be extended to a 
broad range of 2-D gas - 2-D solid phase transitions. In 
addition this type of measurements can give important 
information about dynamics of first order phase transi- 
tions on surfaces and about 2-D interface motion. 
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